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Abstract 

The conventional parametrizations of lattice static interquark force are shown to produce 
a mismatch of coupling constant near the matching point R m ~ 0.2fm. The running coupling 
constant of the background perturbation theory yields instead a selfconsistent description of 
lattice data over all interval of distances, 0.05fm< R <l.l.fm, and is an analytic function with 
the property of asymptotic freedom at small R and it is freezing at large R. 
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1. Recently the interaction betwen static quark and antiquark was intensively studied in 
lattice version of pure gauge SU(3) and SU{2) theories [1-6]. These investigations seem to 
be very important because just in lattice approach both perturbative (P) and nonperturbative 
(NP) effects and their interference can be measured and thus we have an opportunity to get 
an information out of the framework of pure P theory. To this end static potential V(R), or 
force F(R) = 4^, has to be measured over large interval of interquark separations R which 
is possible only on large volume lattices. Such lattice data were presented in ]4], |5| for 577(3) 
gauge theory and here we shall mostly discuss results obtained in [|] (lattice 36 4 at (3 — 6.5, 
a -1 = 4, 13 GeV or lattice spacing a = 0.048 fm) where the force was measured on the large 
interval a < R < 23a or 0.05 fm < R < 1, 1 fm. 

At small distances, R <0.2 fm, the regime of asymptotic freedom was observed and lattice 
data were well described by only P contribution ||. U 

F p (R) = I (i) 

with the effective running coupling constant (c.c.) ap(R) defined by two-loop expression which 
was taken in |4], [j| in the form: 



47T 

aF{R) = m 
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Here b = 11 and b\ = 102 are the usual coefficients for the /3-function with the number of 
flavours rif = 0. 

The "experimental" lattice points a(R) were defined in following way 

where V C (R) is the corrected static potential which takes into account the lattice artefacts 
responsible for the lack of rotational invariance. If we accept that NP contribution is small at 
small separations, R < 4a, then a(R) (3) has to coincide with ap(R) in (1). Below we give 
values of a(R) at two points R/a which are important for our further discussion (the numbers 
are taken from Table 3 in paper jlj): 

R 

- = 3.3839 or R = 0.16 fm, a(R) = 0.317 (10) (2) 
a 

- = 3.9241 or R = 0.19/m, a(R) = 0.333 (6) (1) (4) 

a 
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From the lattice analysis [Q it is clear that these two points R/a lie just near the boundary of 
the region (denote it as R m = 4a = 0.191 fm) where P contribution (1) dominates. At R > 4a, 
when both P and NP contributions are becoming important, another parametrization is widely 
used [4-6], it corresponds to the Cornell potential type of interaction. This force can be written 
(in lattice units) as 

and from lattice fit the following values of K and E were obtained ~ = fjj 

E = 0.278 (7), Ka 2 = 0.0114 (2) (6) 

The string tension K in (5) can be found from (6) and for a -1 = 4, 13GeV it gives \[~K = 
440MeV. Also from (6) it follows that the asymptotic value of c.c. is 

3 

a F (asym) = -E = 0.2085 (5), (7) 

which we can compare with c.c. values (4) near the matching point R m = 4a. From this 
comparison one can see that the asymptotic value (7) is about 40% less than c.c. a(R m ) = 0.33 
which was directly measured on the lattice. Thus if the parametrization (5) is used at R > 4a 
and the purely P force (1) is taken at R < 4a then the resulting function (Xf(R) is discontinuous. 

To overcome this difficulty and get selfcosistent description of static lattice potential or 
force we suggest to use another approach which takes into account the behaviour of running 
c.c. in background vacuum fields. As it was shown by Yu.Simonov in |7|, || this running c.c, 
o>b{R), is an analytical function at all distances R and its expression simplifies at small and 
large distances. At small distances, R~ 2 3> A^, (Xb{R) coincides with standard P coupling 
constant, i.e. manifests the property of asymptotic freedom. At large values of t B (see formula 
(9)) the following expression for olb{R) was obtained, 



a B {R) = t— 1 - ^ — — ), Rm B ^>l (8) 

Here 

, 1 + R 2 m 2 R , , 

and m B is a characteristic background mass which defines the asymptotic (exponential) falling 
off of self-energy part of two-point function Y\(x,y). This mass can depend on the process or 
channel considered and, in particular, in e + e~ annihilation tub coincides with the lowest hybrid 
mass, M/i(0 ++ )||. The lattice measurements 0] has given M/ l (0 ++ ) = 1.5GeV. 

In this letter we analyze a B (R) as obtained from static QQ force. On the lattice the static 
interaction was measured for Eucledean time Tg which is usually not large, e.g. in [|J lattice 
measurements were done at Tg = 3a -j- 5a ~ 0.2fm and for such Tg quark and antiquark are 
situated rather close to each other. In this case background mass m B can be smaller Mh but 
larger than two constituent gluon masses M{2g) = 2Mg rs l.OGeV |§. In our analysis we shall 
use two different values of mg: 

m g) = M h (0 ++ ) = 1.5GeV, m (2) M(2g) = l.OGeV (10) 
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From (8) one can see that this expression, derived at large t B , has a correct P behaviour at 
small R, i.e. has the property of asymptotic freedom. The important second feature of 0£b(R) 
is that it has no infrared singularity and can be applied at any distances. Also the expression 
(8) gives the correct limit at large R when a B (R) tends to a constant value as it was shown in 
0. For this reasons we can use (8) as a convenient interpolation formula over all interval of 
interquark separations R. 

Now to coordinate our analysis with lattice approach we introduce, instead of (8), slightly 
modified expression for o>b(R), 



<*b(R) = — (1 + 15-^1 > (11) 



4-7T / b\ In ts 
b t B V b o f B 

which is going over into the P expression (2) at ms = 0. In (11) i# is defined by (9). From (11) 
one can easily see that olb{R) is freezing at large distances, as(R) — ► const because t B — > const 
if m 2 B R 2 > 1: 



2 

t B CR->oo)— ►t 00 = Zn^- (12) 



R 



The asymptotic value of og(-R) is given by 




a B {asym) = — — 1 + -r~ , °° (13) 



To find a B (asym) we should know, besides ms(10), also Ar. Here we shall choose Ar to be 
very close to those values A R which were obtained in [|J, Ajjg = 256(20)MeV, and in Aj^ 
= 293 (18) MeV. As it is well known |A fi = 1.048A W i.e. A H defining the force and 

Ajg practically coincide within lattice errors. 

Then in the presence of background fields the force F(R) can be presented as a sum 

F(R) = F {PB \R) + F NP (R) } (14a) 
where instead of pure P term (1) we introduce the modified term, 

pPB{R) = 4 a^fl) (M6) 

with olb{R) taken in the form (11). 

Some remarks are needed about our choice of NP part of interaction. We shall consider 
two possibilities: case A . when NP potential is defined by linear potential at all distances or 
F NP (R) = Ka 2 = const as in (5), and case B . when NP interaction is defined through bilocal 
correlators D(x) with a finite correlation length T g pT| . The first measurements of correlators 
and T g were done in [T2| where the value T g = 0.2fm (or 5 = T^ 1 = lGeV) was found. Note 
that the case A corresponds to vanishing correlation length T g , or 5 ^> lGev. The potential 
and force, corresponding to a finite T g and the exponential form of the correlator D(x) were 
calculated in |T3] ], 

F NP (R) = a 2 ^- = 2a 2 J dv [ dX D(V\ 2 + v 2 ) = -Ka 2 (** tKAtUt (15) 
dR J J ix Jo 
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where K\(t) is the modified Bessel function. From (15) one can see that (i) if T g is finite then 
the force is diminishing (F(R) ~ R) at R — > 0; (ii) F(R) is approaching the constant value, 
Ka 2 , only at distances Rf or m ^ 3T 9 ~ 0.6fm. This second statement is in agreement with Bali, 
Schilling result |b| that the formation of string takes place at distances Rf orm ~ 0.6fm. 



Now we give two sets of parameters used for our calculations, 
case A : 

A A = 2A0MeV, m B = 1.5GeV, 
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case B: 



a B {asym) = 0.240; e A = ^a B {asym) = 0.320; F NP (R) = Ka 2 = 0.0114 (16) 



A B = 280MeV, m B = l.OGeV 
a B (asym) = 0.343; e B = -a B (asym) = 0.457; F NP (R) is given by (15) (17) 



3 



with 5 = lGeV or T g = 0.2fm 

Our calculations of the force (the cases A and B) at R < 4a are presented in table 1 together 
with lattice data from M and the P calculations of F P (R) with A R =289 MeV, which slightly 
differs from the value A-^=r = 256 (20) MeV or A^ =268MeV which was predicted in We 
have used A R =289 MeV because for this value we get very good X 2 { N ) f — ^p) whereas for 

A R = 268MeV N * - = is too large. At R < 4a our calculations of F(R) give also a good 

description of lattice data: in case A we get — = ^ and in case B — = %^ is better (see 
table 1). 

At medium and large separations, R > 4a, the values of F(R) (the cases A and B) are 
presented in table 2 together with 15 values found on the lattice. In both cases we have 
obtained a good agreement with lattice data. In case A we get -rj^ — = ^ir and in case B - 

00 ° Nd.o.f. 15 

— * — — -j|P which should be compared with — = obtained in [Q for the parametrization 



M - — = which should be compared with 
(5) corresponding Cornell potential. 

Thus we can conclude that with running c.c. of background perturbation theory it is possible 
to provide a good description of lattice static interaction over all interval of separations R, 
a < R < 23a or 0.05fm< R < 1.1. fm. It holds in both cases, A and B. The c.c. a B (R) 
strongly depends on the parameters m B and A^ (mostly on (t^-) 2 ) but still we are not able 
to distinguish between two sets of these parameters. There are two reasons for it. First, the 
lattice errors are still too large in [Q, jj| and, secondly, the fitted values of m B and A R depend 
indirectly on NP contribution to the force. At present this F NP (R) is not known unambiguously. 
The constant force, F NP (R) = Ka 2 , which is widely used in lattice analysis, corresponds to 
vanishing T g (T g = 0) whereas the finite value of T g = 0.2fm was obtained in ||12|| . So it would 
be important to determine F NP (R) from the lattice measurements in unambiguous way. Then, 
having more precise lattice data on the force and knowing T g with good accuracy, one has an 
opportunity to determine unambiguously m B and A R . 

This a B (R) is an analytical function of R at any distances and freezing at large separations. 
The main difference between cases A and B considered here is in the asymptotic value of 
a B (R) or in constant e: = 0.32 and e B = 0.46, i.e. e B is about 40% larger than e^. 
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It is interesting to note that our is very close to the fitted values of e = 0.31 -j- 0.315 
obtained in || for static potential with linear confining term which corresponds to our case A. 
On the other hand our value of = 0.46 is close to Coulomb parameter e p h en , usually used 
in phenomenological potential models, e.g. to describe heavy quarkonia spectra the values of 
e P hen( n f = 3) ~ 0.45 -T- 0.55 are needed. Notice also that the fitted values of and e# found 
here are about (30 4- 60)% larger than the value of e(string)= ^ = 0.262 predicted in bosonic 
string theory at large distances [[|. 

Author is grateful to Prof. Yu.A.Simonov for fruitful discussions. This work was partly 
supported by RFFI-DFG grant 96-02-00088G. 
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Table 1. Lattice data and theoretical calculations of the force F(R) 
at small distances, R < 4a. 



R/a 


PT 

A R = 289MeV") 


Case A 
A R = 240MeV b ) 


Case B 
A R = 280MeV c ) 


T if\ ft ir*p H p\ t f\ 
from 


1.2071 


0.1662 


0.1623 


0.1639 


0.1607 (55) 


1.7071 


0.0956 


0.0944 


0.0940 


0.0930(23) 


2.1180 


0.0686 


0.0685 


0.0674 


0.0664(48) 


2.5322 


0.0527 


0.0530 


0.0517 


0.0523(6) 


2.9142 


0.0432 


0.0438 


0.0424 


0.0424(57) 


3.0811 


0.0400 


0.0404 


0.0393 


0.0368(49) 


3.3839 


0.0353 


0.0361 


0.0348 


0.0371(14) 


3.9241 


0.0292 


0.0302 


0.0287 


0.0290(6) 



a) calculations in perturbative theory with A R = 289Mev; 

b) m B = 1.5GeV, Ka 2 = 0.0114; 

c) m B = l.OGeV and F NP (R) is given by (15) with Ka 2 = 0.0114, S = lGeV; 

d) lattice errors are found from the paper ||, Table 3 (/? = 6.5; a" 1 = 4, 13GeV). 
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Table 2. Lattice data and theoretical calculations of 
force F(R) at large distances, R > 4a. 



R/a 


Case A a ) 
A R = 240MeV 


Case B b > 
A R = 280MeV 


Lattice data c ^ 
from || 


4.2361 


0.0278 


0.0268 


0.0286(8) 


5.0645 


0.0231 


0.0218 


0.0226(5) 


5.8284 


0.0204 


0.0191 


0.0196(19) 


6.1623 


0.0194 


0.0183 


0.0190(16) 


6.7678 


0.0181 


0.0171 


0.0174(9) 


7.6056 


0.0167 


0.0159 


0.0161(10) 


8.2426 


0.0160 


0.0152 


0.0161(11) 


9.00 


0.0152 


0.0146 


0.0149(3) 


11.00 


0.0139 


0.0136 


0.0130(4) 


13.00 


0.0133 


0.0131 


0.0134(5) 


15.00 


0.0128 


0.0127 


0.0121(6) 


17.00 


0.0125 


0.0126 


0.0132(6) 


19.00 


0.0123 


0.0124 


0.0117(6) 


21.00 


0.0121 


0.123 


0.0125(8) 


23.00 


0.0120 


0.0122 


0.0126(7) 



a) see footnote (b) in Table 1; 

b) see footnote (c) in Table 1; 

c) lattice data for are taken from [|J, Table 1. 
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